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ABSTRACT 

A quadratic  spline  approximation  in  C^  to  the  solution  of  the  Abel 
integral  equation 

^ 1 

j -■  — y(s)  ds  = f (x)  , X 0 , 

0/2  2 

/x  -s 

is  constructed.  It  is  shown  that  if  y'"  is  Lipschitz  continuous,  then  the 
approximation  and  its  first  two  derivatives  converge  to  the  corresponding 
exact  solutions  at  each  point  in  the  interval  of  integration,  the  orders  of 
convergence  being  two,  two  and  one,  respectively.  In  addition,  if 

4 

y £ C [0,1],  then  an  asymptotic  error  estimate  for  the  derivative  of  the 

approximate  solution  is  obtained,  and  convergence  of  the  approximate  solution 

5/2 

is  proved  to  be  of  order  h . The  method  is  illustrated  by  a numerical 
example . 


AMS  (MOS)  Subject  Classif ications : 45EJ0,  45L]0,  65R05. 

Key  Words:  Abel  integral  eqtjation,  Global  approximate  soliJticn,  Quadratic 
splines.  Higher  order  methods.  Asymptotic  error  estimates. 
Product  integration. 
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SIGNIFICANCE  AND  EXPLANATION 


The  Significance  of  the  Abel  integral  equations  in  physical  applications 
and  some  description  of  numerical  methods  of  solution  have  been  discussed  in 
the  Significance  and  Explanation  for  a recent  MRC  report  I#1904]  • In  that 
report,  an  analysis  of  a linear  spline  method  of  solution  was  given. 

In  the  present  paper,  a global  approximation  method  is  developed  for 
obtaining  higher  accuracy  results  involving  higher  order  splines  with  full 
continuity.  The  quadratic  spline  case  is  investigated  in  detail.  The  tech- 


nique is  to  differentiate  the  original  equation,  and  solve  the  differentiated 


equation  by  using  a quadratic  spline  in  c|^.  The  computational  effort 
required  is  only  marginally  greater  than  that  required  for  a linear  spline 
solution  of  the  original  equation.  Convergence  is  obtained  not  only  for  the 
approximate  solution  but  also  for  its  first  two  derivatives  at  each  point  in 
the  interval  of  integration. 


■ 
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A HIGHER  ORDER  GLOBAL  APPROXIMATION  METHOD 
FOR  SOLVING  AN  ABEL  INTEGRAL  EQUATION 
BY  QUADRATIC  SPLINES 


Hinq-Sum  Hung 


1.  Introduction.  A direct  numerical  method  using  quadratic  spline  of  class  C will  be  pre- 
sented and  analyzed  for  solving  the  Abel  integral  equation: 

X 

(1.1)  J y(s)  ds  = f (x)  , X ^ 0 

0/2  2 
/x  -s 

This  is  a clf'ssical  equation,  and  it  occurs  very  often  in  applications  (see,  for  exampxe. 

Noble  [13]  and  Anderssen  [1]).  It  is  well  known  that  this  equation  can  be  converted  to  give 


, , 2 d f sf(s)  , 

= ¥d^  J 

0/2  2 

/x  -S 


= - {f(0)  + X / ds>  , X > 0 . 


0/2  2 
/x  -s 

Methods  based  on  the  numerical  evaluation  of  this  inversion  formula  have  been  considered  in 
[6]  and  [11] . They  all  have  to  deal  with  the  presence  of  a derivative  in  the  inversion  for- 
mula, and  are  inapplicable  to  equations  for  which  an  inversion  formula  is  not  known.  Direct 
methods  for  solving  equation  (1.1)  and  related  equations  have  been  proposed  and  studied  by  a 
number  of  authors  [1-5,  7,  8,  11,  12,  14,  15). 

Recently,  step  by  step  methods  based  on  the  concept  of  product  integration  developed  by 
Young  [16]  for  the  solution  of  (1,1)  by  using  linear  splines  have  been  considered  and  analyzed 

by  several  authors  ([2],  [3],  [8]).  Atkinson  [2]  gives  a convergence  theorem  but  does  not 

2 2 

prove  that  the  convergence  is  0 (h  ).  Benson  [3]  obtains  0 (h  ) and  also  derives  an  asymp- 
totic formula  for  the  discretization  error,  but  his  method  depends  on  a comp 1 i oa t ed  analysis 
of  product  integration.  Hung  [8]  obtains  similar  results  as  Benson  |3I  uncer  slightly  weaker 
conditions  and  by  a much  simpler  method. 
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It  might  be  supposed  that  the  use  of  splines  of  order  greater  than  one  would  result  in  a 


higher  order  approximation  method  for  equation  (1.1).  Such  an  approach,  however,  fails, 
because  of  instability,  in  the  case  that  the  spline  is  of  full  continuity  and  the  approxima- 
tion to  the  solution  is  based  on  using  equation  (1.1)  alone  (see  Appendix).  The  divergence  of 
higher  order  spline  methods  arise  from  the  strictness  of  the  continuity  requirements  on  the 
spline  functions.  It  is  possible  to  generate  higher  order  methods  by  relaxing  the  continuity 
i requirements  on  the  spline  functions,  and  this  is  the  type  of  methods  considered  by  Brunner 

and  Linz.  Brunner  [4,  5]  uses  piecewise  polynomials  of  a given  decree  and  of  class  C to  gen- 
erate approximate  solutions  for  equation  (1.1)  as  well  as  equations  of  a similar  type,  and 
shows  that  convergence  of  order  m is  obtained  if  m*"^  degree  piecewise  polynomials  are 
used.  His  method  is  based  on  using  the  original  equation  alone,  and  is  basically  a bloc)c  by 
bloc)(  method  in  the  sense  of  (10).  Linz  [11]  shows  numerically  that  a fourth  order  Hermite 
interpolation  method  can  be  constructed  for  equation  (1.1)  by  using  cubic  spline  in  . His 
I method  is  based  on  using  the  differentiated  form  of  equation  (1.1)  simultaneously  with  equa- 

! tion  (1.1) . 

In  this  paper  we  show  that  the  divergence  of  higher  order  splines  with  full  continuity 
depends  upon  the  particular  method  used,  and  is  not  necessarily  a result  of  strictness  of  t)ie 
continuity  requirements.  We  develop  a global  approximate  solution  for  equation  (1.1)  by 
quadratic  splines  in  . The  method  to  be  described  in  Section  2 is  self-starting  and  is 
based  on  using  the  differentiated  form  of  (1.1)  alone.  In  Section  3,  wo  prove  convergence  for 

the  approximate  solution  as  well  as  its  first  two  derivatives.  A simple  asymptotic  error 

I 

I formula  is  derived  in  Section  4 for  the  derivative  of  the  approximate  solution  by  which  con- 

I 5/2 

vergence  of  the  approximate  solution  itself  is  shown  to  be  of  order  h . Finally,  a numeri- 
cal example  is  presented  in  Section  5. 

In  this  paper,  it  is  assumed  that  0 < x < 1,  but  this  restriction  is  not  essential. 


I 

) 


I 


2.  The  (Juadratic  Spline  Method, 
constant  step  size.  Let  y(x! 
approximations  to  y(x.l,  y'(x.) 


Let  x^  = ih,  i - O,  1,  2,  •*•,  where  h is  an  arbitrary 
be  the  exact  solution  of  (I.l),  and  let  Y.,  Y’  denote 


1 

respectively.  We  use  a quadratic  spline 


P{x)  in  C^[0,11 


1 


T 


-2- 


with  knots  at  the  points  as  an  approximation  to  y{x),  i.e.,  for  i = 0,  1,  2, 

(2.1)  P(x)  = + I i(2u^(x)  - u^(x))y^  + u^(x)Y^^^l  , 1 X 1 . 

where  (x)  = (x-x^)/h.  Differentiating  (2.1)  we  obtain  the  derivative  of  P(x), 


(2.2)  P'(x)  = u.(x)Y’.  , - u.^,  (x)Y',  , X.  < X < x_,  . 

1 1+1  1+1  1 1 “ 1+1 

Both  P(x)  and  P'(x)  are  continuous  at  the  knots. 

In  order  to  obtain  a numerical  solution  of  (1.1)  by  using  the  quadratic  spline  (2.1).  we 

have  to  introduce  the  differentiated  form  of  (1.1).  Assume  that  we  can  differentiate  (1.1). 


(2.3)  / - y'  (s)  ds  - xf  ’ (X)  , X > 0 . 

0 Jl  2 

/x  -s 

The  approximate  solution  P(x)  of  the  integral  equation  (1.1)  is  derived  from  values 

P(x.)  = Y.  and  P'(x.)  = Y*.  . We  know  how  to  find  Y*  , assuming  that  we  know 
i 1 i i K+i 

Y| , i = 0,  1 , • • • , k.  We  can  then  deduce  from  the  equation  obtcd.ned  by  setting 

X = in  (2.1),  which  gives 

<2-4)  Vi  = \ " I " "kn>  • 

To  compute  Y^,  we  require  that  P*  (x)  satisfies  (2.3),  i.e.,  y' (x)  is  replaced  by  P' (x) 
derived  from  the  values  P'(x^)  “ Y'^  computed  from 


0/2  2 
/x,  -s 
k 


P'  (s)  ds  = Xj^f ' (Xj^) 


This  can  be  rewritten  in  the  form: 


I “k.i’^i ' 

i»0 


k = 1,  2, 


i 

! 


L.ik: 


k,0 


-/ 


x_  / 2 

0 /x,  -s 
k 


Uj^{s)  ds 


= -^U..^(S)  ds  - / 


i+1 


X.  , 2 

X-l/Xj^-S 


/2  2 
/x,  -s 
k 


i=l,---,k'l  , 


\,k  = / ^ 


X , /T~  2 
k-l/x,  -s 
k 


u,  , (s)  ds  . 
k-1 


Equation  (2.6)  is  a nonsingular  triangular  system  for  Y\  because 


(2.8) 


for  k = 1,  2, 


• 4 1 5k-l 


r 4 ■,  bk-i  , C 4 


Since  y(0)  = — f(0)  and  y'(0)  = f'(0),  we  take  = y(0)  and  Y^  = y'(0).  Tfie 
values  Y^,  Y^  (i  = 1,  2,  •••)  can  tlien  be  determined  successively  by  (2.6)  together  with 

(2.4).  Note  that  this  approach  requires  no  starting  procedure. 

An  estimate  of  y"(x)  is  given  by  the  second  derivative  of  (2.1).  If  we  denote  this 

(constant)  estimate  of  y"  (x)  in  x^  £ x < x^^^^  by  YV,  this  gives 


(2.9) 


Y"  = r i'i-  , ' Y!) 
1 h 1+1  1 


X . < X < X 


i+1 


3.  Convergence  Pesults.  The  proofs  of  Theorem  3.1  and  Theorem  4.1  require  the  following 
lemma  which  is  a consequence  of  the  standard  results  for  regular  infinite  systems  of  algebraic 
equations  by  Kantorovich  and  Krylov  in  [7,  p.271 . 

Lemma  3.1.  If  there  exist  a constant  c ^ 0 and  an  integer  N > 1,  all  independent  of  k, 


such  that 


|xj  < C 


k + l ' 


i = 0,  1 , • • •,  N 

k 

- ^ J Is  I , k = N,  N + 1,  ••• 

i=()  ^ 


with 


k = N,  N + 1,  • - • 


1 - I |a 


i=0 


k+l,i' 


> 0 


then 


^«k 


|x. 1 < C 


0,  1,  2, 


The  proof  of  this  lemma  can  be  found  in  [8,  p.3] . 

Let  y(x)  be  the  exact  solution  of  (1.1)  and  define  the  discretization  error  function 
€(x)  by  £ (x)  = y(x)  - p(x) , where  p(x)  is  the  quadratic  spline  approximation  to  y(x) 
obtained  from  our  numerical  method.  Denote  £ by  £|*^'  (r  = 0,  1,  2).  We  state  the 

following  theorem: 

Theorem  3.1.  If  y"'(x)  is  Lipschitz  continuous  on  [0,1],  then 


(3.1)  |£|^|  = 0(h  ) , k = 1,  2,  ••• 

Proof ; Applying  Taylor's  theorem,  it  is  not  difficult  to  show 


(3.2)  y(x)=y(x^)+|  [ (2u^  (x)  - u^  (x)  )y ' (x^)  + u?  (x)y ' (x^^^)  ] + (x)  , 1 ^ 1 


where 

(3.3) 


'^(x)  = — 


*i+l 


/ y"'(s)  ds  + u^(x)  / u^^^  (s)y (s)  ds 


with  u^ (x)  as  defined  in  (2.1). 

Subtracting  (2.1)  from  (3.2)  gives 

(3.4)  £(X)  = £ + ^ [(2u  (X)  - u^(x))£!  + u^(x)£!^,)  + ^(x) 

Differentiating  (3.4),  we  have  the  derivative  of  £ (x) , 

(3.5)  £'(x)  = [u.(x)£'  -u.  ,(x)£!1  +v>'(x) 

where  ifi ' (x)  can  be  represented  in  Lagrange  form  as 

^2 


X.  < X < X.  . , 
1 — — 1+1 


X.  < X < X . . , 
1 — - 1+1 


(3.6) 


■^Mx)  = — y"’ (n^(x))u^(x)u^^j^(x) 


X.  < n. (x)  < X.  , 
1—1  - 1+1 


Since  both  y' (x) 

and 

P'  (x) 

satisfy  (2.3) 

at  each  knot  x = x, 
k 

k-1 

’'i+l 

(3.7)  1 

/ 

s 

£'  (s)  ds  = 0 

, k = 1,  2,  ••• 

i=0 

X . 

/l  2 

1 

/x,  -s 

V 

-5- 


This  can  be  rewritten  as 


i=0 


, , 5C  . , 

k-1  1+1 


\ = / 


k = 1,  2, 


ip ' (s)  ds  , 


i=0  X.  2 2 

1 vx,  -s 
k 

and  the  w,  . 's  are  defined  in  (2.7). 
k,i 

Multiply  (3.8)  by  k,  difference  the  resulting  equation  for  k and  k + 1,  and  then 
divide  by  (k  + k+1  required  error  equation 


>+1  \+l,i®i 


kWk  i - 


(k+l)w, 


k = 1,  2, 


i — 0/  X,  ***r  r 


k+1, k+1 


(k+l)R^^l  - kR^ 


^k  (k+l)w 


k+1, k+1 


k = 1,  2, 


Note  that  (3.10)  is  not  the  only  possible  error  equation  that  can  be  constructed  for  the  proof 
of  convergence.  For  instance,  we  can  limply  difference  equation  (3.8)  for  k and  k+1. 
However,  the  procedure  used  here  simplifies  the  asymptotic  error  analysis  in  Section  4. 
Equation  (3.10)  implies  that 


i");+ii  ^ Jq  i^+i,ii  ^ 


k = 1,  2, 


Since  ~ it  is  easily  shown  from  equation  (3.8)  by  using  Lemma  3.3(a)  below  that 
2 

£|  = 0(h  ) for  i = 1,  •••  , K,  with  K as  defined  in  Lemma  3.2  below.  On  the  basis  of 
this  together  with  Lemma  3.2(b)  and  Lemma  3.3(c)  below  we  can  apply  Lemma  3.1  on  (3.13)  and 
obtain  (3.1).  Hence  the  proof  of  Theorem  3.1  is  completed. 

Theorem  3.2.  Let  the  assumption  of  Theorem  3.1  be  satisfied,  then  for  k = 1,  2, 

(3.14)  = O(h^)  . 

Proof : Setting  x = ’'i+i  (3.4)  and  noting  that  + = O(h^),  we  obtain 


£.  =«.  + ^ (£:+£!,)  + O(h^)  . 

i+1  1 2 1 1+1 


I i 


I f 


Since  = 0,  by  means  of  (3.1),  (3.14)  inmediately  follows  from  (3.15)  by  induction. 

In  the  following  theorem  we  can  show  that  the  approximate  solution  P(x)  along  with  its 
first  and  second  derivatives  converges  to  the  corresponding  exact  solutions  at  each  point  in 
the  interval  of  integration  [0,1] . 

Theorem  3. 3.  Let  the  assumption  of  Theorem  3.1  be  satisfied,  then  for  any  fixed  x e [0,1], 

|e(x)|  = 0(h^)  , |£*  (x)  I » O(h^)  , | e"  (x)  | = 0 (h)  . 

Proof;  By  means  of  (3.1)  and  (3.14),  Theorem  3.3  follows  immediately  from  equation  (3.4), 
(3.5)  and  the  equation  resulting  from  differentiating  (3.4)  twice. 

Lemma  3.2.  If  the  defined  by  (3.11),  then  there  exists  an  integer  )c  ^ 1, 

independent  of  h amd  )c,  such  that 

(a)  10  , i-0,  l,-",)c-l,k>,l, 

\tl,lc 

k 


Ic  > K 


(b) 


Proof  of  (a) ; From  (3.11),  using  (2.7)  and  (2.8)  we  have 


®)c+l,0  ” ()t+l)w, 


/ 


lc+l,k+l  X 


A+1  ■ 


s (Xj^-s) 


ds  > 0 


since  the  integremd  is  non-negative . Similarly,  we  cem  prove  that  i — ° i = 1, 

•••,  k - 1.  Finally,  it  is  easy  to  show  that  tends  to  3 - 2^2  as  k increases, 

therefore  there  exists  an  integer  K,  independent  of  h and  k such  that  k — ° 

k ^ K. 

Proof  of  (b)  ; Since  by  using  (2.7)  we  obtain 

k «k 

,L"k,i-(,  -7='“-"k 


(3.16) 


i»0 


0/2  2 
/Xk-s 


for  k = 1,  2, 
(3.17) 


',  it  cem  easily  be  shown  that  for  k = 1,  2,  ••• 
k 

.^„\+l,i  ” ^ (k+l)w, 


(2k-H)h 


i=0 


k+l,k+l 


-7- 


si-M 


k+l 


be,... 


-r'.'Wi 


Let  Lj  be  the  Lipschitz  constcUit  for  y'"  . Ttien  by  straightforward  estimation  aiid  noting 
that  hk  ^ 1,  we  obtain  from  (3-19) 

<3-20)  IVl  - ' 

where  (L^  + 3Mj) . 

Proof  of  (c) ! From  (3.12),  using  (3.18),  (3.20)  and  (2.8),  it  can  easily  be  shown  that  for 


k = 1,  2,  •••, 


^\+l  " \ ^ k+1  ^ 


l\l  1 ^ 


k+l,k+l 


-^3  ^ ’ 


where  C3=-(C^+C2). 


4.  An  Asymptotic  Error  Fomula.  In  this  section,  we  obtain,  under  the  cissumption  that 

4 

y e C [0,1],  an  asymptotic  error  eiqpression  for  which  is  essential  in  proving  that  con- 

5/2 

vergence  of  the  approximate  solution  is  of  order  h 
4 

Theorem  4.1.  If  y £ C [0,1],  then  for  k 1,  2,  ••• 

h2  2 

(4.1)  £j^  = ^ y’"  (Xj^)  + 0(h  //It)  . 

4 

Proof ; Since  by  assun^tion  y (x)  € C [0,1]#  it  is  not  difficult  to  show  that,  for 
X £ [X.  , X.^3] 

(4.2)  £(x)  = + 7 [(2u^(x)  - u^(x))£^  + V (*i  > t2u^  (x)  - 3u^  (x)  ] + p (x)  , 


p (x)  = ^ [2  J Y^^'(s)  ds  - 3u^(x)  / u^^j^(s)y  (s)  ds] 

X.  X. 

1 1 

with  u.  (x)  as  defined  in  (2.1). 

1 


Differentiating  (4.2)  we  obtain  the  derivative  of  e(x), 


(4.3)  £*  (x)  = [u^  (x)£^^3  - u^^3  (x)£|]  +— y'"  (x^)u^  (x)u^^j^(x)  +p' (x) 

where  p ’ (x)  can  be  represented  in  Lagrange  form  as 


1 1 < 


(4.4)  p'  (x)  = ^ (n(x))  - (x)y  (n  ) ) , x^  f.n(x)  f.  x^^^^  . 

.2 

Adding  and  subtracting  - — [u^  (x)y'"  (x^^^^)  - u^^j^  (x)y'"  (x^)  ] to  the  right-hand  side 
of  (4.3)  and  using  the  Taylor  series  expansion  for  y'"  (x)  at  x^  for  x = 
difficult  to  show  that 

^2 

(4.5)  £'  (x)  = [u^(x)  y- (x^^j^))  -u^^^(x)  (e[ -^y"'  (x^))]  +$(x)  , x^^f^x^x^^^ 

where 


2 3 

♦(X)  y'"  (x^)  [u^  + [p'  (x)  Uj^  (x)y  (C(x^^j^))  ] , ^5  (x)  £x. 


i+l 


By  substituting  (4.5)  into  (3.7)  we  obtain 


(4.6) 
»diere 

(4.7) 


J^''k,i<^i  - 12  y'"(*i))  = \ 


k = 1,  2, 


lc-1  *i+l 


t (s)  ds 


i=0  X. 
1 


J 2 2 

/x^-s 


Comparing  (4.6)  and  (3.8)  we  note  that  Rj^  is  converging  faster  than  If  we  derive 

an  error  equation  in  the  same  way  as  in  the  proof  of  convergence  in  Section  3 without  modify- 
ing (4.6)  we  expect  the  parts  in  the  derived  error  equation  vdiich  correspond  to  Sj^  and  pj^ 
in  Lemma  3.1  would  have  unbalanced  rates  of  convergence,  with  the  former  converging  faster 

than  the  latter.  By  applying  Lemna  3.1  to  such  an  equation  we  will  fall  to  obtain  the  result 

h^ 

we  expect..  To  avoid  this,  we  define  = /i+l  (c^  - jy  y'"(x^))  and  rewrite  (4.6)  as 


(4.8) 


1,  2, 


iio  ^ 

Now  multiply  (4.8)  by  k,  difference  the  resulting  equation  for  k and  k + 1,  then 
divide  by  (k  + l)w.  /i(k+2  to  yield  the  required  error  equation: 


(4.9) 


^k+1  = A+i,i^i  ^ K 

1=0 


k = 1,  2, 


where 


-10- 


®lc+l,)c  - ° ' 

(b)  1 - I -7=  ' 

i=0  .W 


i = 0,  1,  • • • , X - 1,  K > 1 


X > K , 


k > K , 


for  em  appropriate  C,  0 ^ C < 1. 

Proof  of  (a):  By  means  of  Lemma  3.2(a),  part  (a)  follows  iamediately  from  (4.10). 

k 

Proof  of  (b) ! Since  it  would  be  very  coii(>lex  if  we  estimate  ^ ®X+1  i intro- 


Vi.i= 


(4.15)  “ *k-fl,i  (k-H)w, 


0 < C < 1,  1 “ 0,  1,  •••,  k 


X-tl  ,k-*-l 


By  using  Lemna  3.2(a)  and  noting  the  non-negativity  of  the  Wj^  ^ from  (2.7),  we  obtain 


(4.16) 


>0  , i - 0,  1,  •••,  X - 1,  X > 1 , 


*k-*-l,k  — 


> 0 , k > K 


where  K is  the  same  K as  defined  in  Lemma  3.2. 

using  (3.16)  euid  (4.16)  it  can  easily  be  verified  t)iat 


Jo  ' 


k > K . 


If  we  can  show  that  for  an 


appropriate  C,  0 < C < 1,  a^^^  . - . 1 


. >0,  then  the  proof 


is  conpleted.  Tto  do  this,  it  is  sufficient  to  show  that  for  i « 0,  1,  •••,  k, 

vi.i - & *"1  •a»k,i '] ['"i-.i - '"♦““k.i.i]  ■ 

By  using  (2<7)  and  letting  s = ht,  it  is  easily  verified  that  for  i * 1,  •••  , k - 2, 


(i-H-t)  dt 


1 -t-  C(1  + i)  / 


\ .2 


since  for  i-l<t<i+l,i*l,  and  0 < C < 1 


ca  + i)  - (i+l)^3  + C’'^  ■ 

- (i-1)^] 


it  is  obvious  fram  (4.18)  that  i — ® i“l.  •••,k-2.  For  simplicity,  t2Uce 

C = 0.  Then  with  this  value  of  C,  it  is  not  difficult  to  show  that  as  k increases 

• (i  - 0,  k - 1,  k)  tends  to  ^h,  (3«^  - 4^  f l)h,  and  ^ (2/2  - Dh. 

K^XfX  OK  3 3 

respectively,  and  are  therefore  greater  than  zero  for  sufficiently  large  k. 

Thus  for  i • 0,  1,  •••,  k,  there  exists  an  integer  K ^ 1,  independent  of  h and  k, 
such  that 

\+l,i^\+l,i  ' • 

By  means  of  peurt  (a),  (4.16)  and  (4.17),  it  follows  that 


^ ” 1^0  ^ ' -^o  “ XIT 

i»0  i“0  >Tc+l 


for  C >•  0,  and  k ^ K,  with  K - max(K,  K). 

4 ... 

Lemma  4.2.  If  y e C 10,1],  then  there  exist  constants  Cj^,  C^,  > 0,  independent  of 

h and  k,  such  that 

(a)  |\|  < •ic  h^  , 

IVl-\li^2^  ' 

(c)  |b^|  < . 

for  k = 1,  2,  where  ^ and  bj^  are  defined  in  (4.7)  and  (4.11),  respectively. 

Proof  of  (a) ; By  repeated  integration  by  parts,  it  is  not  difficult  to  show  that  for 
i = 0,  1,  •••,  k - 1 

.2  *i+l 

(4.19)  y'”(x^)[u^(s)u^^^(8)  t g-1  ds 


X.  /2  2 

r /x^-s 


. Xj  , 2 

h^  f *k= 


(/x?-s^) 


J y'"  (Xj^)u?  (s)u^^j^  (s)  ds 


Using  (4.19)  we  ceui  rewrite  (4.7)  as 


- ;(1)  ^ ;(2) 
\ = \ ^ \ 


]c  - 1,  2, 


,, , .4  )c-l  ^i+1  x?s  , 

*i  [/f?} 


/-»v 

■ - .L ! 


*i+l  ^3  ... 

I ! — = Ip' (s)  + ^ u (s)y'’' (e(x  ))]  ds  , 

i*0  X.  J 2 2 1 i A 


r -'x^- 


with  p'  (s)  as  defined  in  (4.4) . 


Let  Mj  = ^ e”!o,l]  “4  ' X e”^0,ll  straightforward 


estimation,  and  noting  that  h)c  1,  we  obtain  from  (4.20) 


(4.21) 


I I £ M J *ir  h^  + M^jch^ 

^ Cj^  vie  h^  , 


k = 1,  2, 


where  C,  = M,  + M,. 
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Proof  of  (b) ; Subtraction  of  (4.20)  from  (4.20)  with  )c  replaced  by  )c  + 1,  then  by 
straightforward  estimation,  and  noting  that  h)t  ^ 1,  it  is  not  difficult  to  show  that 
for  )i  = 1 , 2 , • • • 


< (13M,  + M.)  + 4M.  h^ 

3 4^4 


< C ^ 

^ ' 


where  “ ^^3  ■*■  ^*^4* 

Proof  of  (c)  : From  (4.11),  using  (4.21),  (4.22),  and  (2.8)  we  obtain  for  )c  = 1,  2,  ••• 


|b  I < ^ 


[ft  * DlS^^i  - 5^1  ns^G  , 


vtfiere  C^  = j(C^  + C^). 


5.  A Numerical  Example.  The  quadratic  spline  Eoethod  described  in  Section  2 was  applied  to 


the  following  Abel  integral  equation: 

X 

/ y <s) 

0 /22^  " 

I'X  -s 

The  exact  solution  is  y (x)  > cos  (x) . 

In  Table  S.l(a)  and  Tedile  5.1(b)  we  list  the  error  e (x) , £' (x)  emd  £"(x)  at  knots  2uid 

at  mid-points  between  the  knots,  respectively  on  (0,31  for  different  step  sizes  h.  The 

5/2  2 

error  £ (x) , £' (x)  amd  £"(x)  satisfy  the  predicted  h , h , and  h dependence,  respec- 

tively. Note  adso  in  Table  5.1(b)  that  the  errors  e“(x)  at  the  mid-points  are  actuadly 
2 

0(h)  although  we  have  not  proved  that  this  will  be  the  case. 

In  TeUsle  5.2  we  list  the  actual  error  (column  2)  of  the  derivative  of  our  approximate 
solution,  i.e.  £'(x),  amd  its  theoretical  estimate  (column  3)  confuted  from  equation  (4.1) 
at  knots  on  [0,3]  for  h = 0.01. 


r » 


6.  Conclusion.  A global  approximation  method  of  order  2.5  to  the  solution  of  the  Abel 
integral  equation  (1.1)  by  using  a quadratic  spline  in  has  been  presented. 

(1)  The  method  is  self-stMting  and  the  step  size  h can  be  changed  at  any  knot  of  the 
quadratic  spline  without  added  coup li cation . 

(2)  It  is  a single  higher-order  method  that  requires  a conpeiratively  small  conputational 
effort.  At  each  step,  only  one  equation  has  to  be  solved,  whereas  other  higher-order  methods 
(see,  for  exaiqple,  [4]  and  [11])  require  solution  of  a system  of  equations  whose  coefficients 
involve  a much  larger  number  of  integrals  to  be  evaluated. 

(3)  Convergence  for  the  approximate  solution  as  well  2is  its  first  two  derivatives  is 
obtained  by  our  method,  the  orders  of  convergence  being  5/2,  2 and  1,  respectively.  In 
addition,  an  asymptotic  error  estimate  is  derived  for  the  derivative  of  the  approximate 
solution. 

(4)  The  method  is  economical  when  the  values  of  the  solution  and  its  derivatives  are 
required  at  a large  number  of  points  where  usual  discrete  methods  of  confutation  will  be  time 
consuming. 

(5)  The  method  should  be  particularly  useful  if  the  derivative  of  f (x)  which  is  on 
the  right-hand  side  of  (1.1)  c2m  be  confuted  analytically  and  is  sinf  le . If  the  derivative  of 
f (x)  c6mnot  be  evaluated  analytically  in  a convenient  way,  or  if  it  is  given  in  tabuleu:  form, 
one  cem  use  appropriate  finite  differences  to  perform  the  approximation.  This,  eis  one  can 
show  easily,  does  not  affect  the  convergence  of  our  method. 

(6)  The  results  developed  in  this  paper  give  some  idea  of  the  necessary  tools  emd  the 
possible  results  for  other  higher-order  spline  approximation  to  the  solution  of  some  more 
general  Abel  integral  equations,  e.g.,  cm  approximate  solution  can  be  obtained  involving  ctibic 
splines  in  C^  by  using  the  equation  resulting  from  dif ferentiati.->q  the  original  equation 


APPENDIX 

Diveraence  of  Methods  Aoolvinq  to  Equation  (1.1) 


Bv  Meems  of  Splines  of  Order  ^ 2 with  Full  Continuit 


For  illustration,  apply  the  method  by  using  a quadratic  spline  in  C 10,1]  to  the 
following  Abel  integral  equation: 

X , , , 


/ • ■■ — y (s)  ds  « 


0/2  2 
*x  -s 

whose  exact  solution  is  simply  y(x)  = x^.  The  divergence  of  the  numericed  results  can  be 
seen  very  clearly  in  Table  A. 

Ted>le  A 

Divergence  of  Numerical  Results 

Bv  the  Quadratic  Spline  Method  Based  on  Using  Equation  (1.1)  alone 


X 

c (x) 

h 

h/3 

e' 

h 

(X) 

h/3 

e"(x) 

h 

1 J 

h/3  1 

0.0 

.OOOOE+0 

.OOOOE+0 

.OOOOE+0 

.OOOOE+0 

-.1698E+0 

.5659E-1 

0.2 

-.  3016  E- 3 

- 

.2168E-3 

-.1208E-1 

-.36  31E-1 

.2  311E+0 

. 3317E+1 

0.4 

-.1338E-2 

* 

.1784E-1 

-.7070E-1 

-.  3044E+1 

. 204  3E+1 

. 2816E4'3 

0.6 

-.585  3E-2 

- 

.1456E+1 

-.3268E+0 

-.2488E+3 

.9952E+1 

.2300E+5 

0.8 

-.2550E-1 

- 

.1184E+3 

-.1443E+1 

-.2024E+5 

.4440E-)-2 

.1871E+7 

1.0 

-.1109E+0 

- 

.9606E+4 

-.6299E+1 

-.1643E+7 

.1942E+3 

.1518E+9 

€ (x)  “ Discretization  error  function. 
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